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Abstract 

The distribution of the azimuthal angle of charged and neutral hadrons relative 
to the lepton plane has been studied for neutral current deep inelastic ep scat- 
tering using an integrated luminosity of 45 pb -1 taken with the ZEUS detector 
at HERA. The measurements were made in the hadronic centre-of-mass system. 
The analysis exploits the energy-flow method, which allows the measurement to 
be made over a larger range of pseudorapidity compared to previous results. The 
dependence of the moments of the azimuthal distributions on the pseudorapidity 
and minimum transverse energy of the final-state hadrons are presented. Al- 
though the predictions from next-to-leading-order QCD describe the data better 
than do the Monte Carlo models incorporating leading-logarithm parton show- 
ers, they still fail to describe the magnitude of the asymmetries. This suggests 
that higher-order calculations may be necessary to describe these data. 
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1 Introduction 



The description of the hadronic final state in deep inelastic scattering (DIS) is influenced 
by perturbative Quantum Chromodynamics (pQCD) in several ways that can be calcu- 
lated through exact matrix elements or leading-logarithm parton showers. Measurements 
of the azimuthal distribution of hadrons in the semi-inclusive process e + p — > e + h + X 
in DIS are sensitive to predictions of pQCD. The azimuthal angle, <fr, is defined in the 
hadronic centre-of-mass (HCM) frame as the angle between the hadron-production plane 
and the lepton-scattering plane as illustrated in Fig. ^ 

Asymmetries in <fi result from the final-state hadron having transverse momentum with 
respect to the colliding virtual photon and the incoming proton. In pQCD, higher-order 
QCD processes such as QCD Compton (QCDC) (7*g — > qg) and boson-gluon fusion 
(BGF) (7*<7 — > qq) are the main sources of these hadrons. These two processes have 
different cj) behaviours [1] as well as a different pseudorapidity, r], dependence, defined 
here with respect to the incoming proton direction in the HCM frame. Figure a) shows 
that hadrons from BGF and QCDC dominate over quark-parton-model (QPM) (7*^ — > q) 
events in the region —4 < ry HCM < 0. In addition, gluons and quarks from the QCDC 
process have different pseudorapidity dependencies, as illustrated in Fig. I2fb). 

The azimuthal dependence for semi-inclusive neutral current (NC) DIS can be written [2- 
5] as: 

da ep ~* ehX 

= A + B cos </> + C cos 2(f) + V sin </> + £ sin 20. (1) 

d(p 

The azimuthal asymmetries, specified by the parameters B, C, T> and S, are extracted 
from the data by calculating the statistical moments of the experimental distributions: 

(cos 6) = — r; (sin. 6) = — -: 

(cos20) = 7^-; (sin 20) = ^j. 

Equation ([!} results from the polarisation of the exchanged virtual photon. The coefficient 
B originates from the interference between the transversely and longitudinally polarised 
components; the coefficient C is due to the interference of amplitudes corresponding to 
the +1 and —1 helicity parts of the transversely polarised exchanged boson. The coeffi- 
cients T> and £ arise from parity-violating weak interactions or longitudinal polarisation 
of the initial lepton beam [3]. They vanish for purely electromagnetic interactions with 
unpolarised beams. 
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It has been proposed [4] to analyse the asymmetry as a function of the transverse mo- 
mentum cutoff, p^, of a detected hadron. Such a cut is efficient in removing QPM 
events. Consequently, at higher p™* values a better agreement should be obtained with 
the perturbative QCD predictions. A model using a resummation formalism [6] to pre- 
dict azimuthal asymmetries has also been proposed. This model predicts that logarithmic 
corrections due to soft parton emission could be large. A recent paper [7] showed that a 
part of the asymmetries previously measured by the ZEUS collaboration [8, 9] may come 
from terms that are not included in the perturbative gluon radiation but are related to 
the intrinsic transverse motion of quarks. Predictions of these models are addressed in 
this paper. 

For NC DIS with an unpolarised lepton beam, the (cos0) and (cos 20) values have been 
measured by the ZEUS collaboration [8, 9] to be at the few percent level. The first 
publication [8] measured the azimuthal distribution for charged hadrons, whereas the sec- 
ond [9] was performed for jets of high transverse energy. The present analysis used the 
energy-flow method, which permits both neutral and charged hadrons to be included in 
the measurements. This analysis was performed using a similar data sample but in an ex- 
tended kinematic range compared to previous publications. In particular, the polar-angle 
range of the measurements was increased with respect to the previous studies [8,9]. The 
energy-flow method enhances the contribution of leading hadrons since the direction of 
each particle in the final state is weighted with its transverse energy [10,11]. Additionally, 
the values (sin 0) and (sin 20) were determined, although they are expected [3, 5] to be 
much smaller than (cos 0) and (cos 20) . 

This paper presents measurements of (cos0), (cos 20), (sin0) and (sin 20) as a function 
of the pseudorapidity, r] ncM , and minimum transverse energy, (instead of pjP [4] ) , 

of the final-state hadrons. The results are compared to theoretical expectations. 

2 Data sample 

The experimental results are based on the data collected in 1995-97 with the ZEUS detec- 
tor at HERA. Protons of 820 GeV collided with 27.5 GeV unpolarised 1 positrons. Neutral 
current DIS events were selected from data corresponding to an integrated luminosity of 
45 pb -1 . A detailed description of the ZEUS detector can be found elsewhere [13]. A 
brief outline of the components that are most relevant for this analysis is given below. 

The high-resolution uranium-scintillator calorimeter (CAL) [14] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 

1 The positrons were transversely polarised due to the Sokolov-Ternov [12] effect, but had no longitudinal 
polarisation. 
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is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a(E)/E = 0.18/ VE for electrons 
and a(E)/E = 0.35/ \[E for hadrons, with E in GeV. 

Charged particles are tracked in the central tracking detector (CTD) [15], which oper- 
ates in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD 
consists of 72 cylindrical drift chamber layers, organised in 9 superlayers covering the 
polar-angle 2 region 15° < 9 < 164°. The transverse-momentum resolution for full-length 
tracks is a{p T )/p T = 0.0058p T © 0.0065 © 0.0014/p T , with p T in GeV. 

The selection criteria were based on an earlier ZEUS investigation [8]. The main require- 
ments on the event were: 

• an identified scattered positron with energy E' e > 10 GeV. Energy deposits in the 
CAL, within an rj — <p cone of radius 1 and consistent with being a photon, were 
removed from the calculation of the positron energy; 

• 100 < Q 2 < 8000 GeV 2 , 0.2 < y < 0.8 and 0.01 < x < 0.1. The quantities x,y and Q 2 
are respectively x-Bjorken, the inelasticity, and the negative square of the exchanged 
boson virtuality. The double angle method was used to reconstruct these variables 
and so determine the direction of the exchanged boson [16]. The mean values of x, y 
and Q 2 were respectively 0.0222, 0.351 and 700 GeV 2 . The final sample consisted of 
16472 events. 

Charged and neutral final-state particles were reconstructed using a combination of track 
and calorimeter information that optimises the resolution of the reconstructed kinematic 
variables [17]. The selected tracks and calorimeter clusters are referred to as energy flow 
objects (EFOs). The EFOs were required to satisfy the following: 

• transverse momentum in the laboratory frame pt > 0.15 GeV; 

• polar angle 9 > 8°. 

These cuts ensured the analysis was performed in a region of high acceptance and high 
detector efficiency An average of about 18 EFOs per event were reconstructed, yielding 
a total of 293 000 objects, each of which provided a value of the azimuthal angle, <j), used 
for further analysis. 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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3 Method 



The energy-flow method was proposed [18,19] to ensure that the inclusive variables such 
as those measured here are infra-red and collinear safe due to the direction of each particle 
being weighted by its transverse energy. 

In this analysis, the moments of trigonometric functions of the azimuthal angle are cal- 
culated as follows. For a function, F(n0 HCM ), where F(n0 HCM ) can be sin(n0 HCM ) or 
cos(n0 HCM ) and n = 1 or 2, the mean value was determined using the formula: 



771HCM 77i /" t ,j.HCM\ 
(f(n ,HCM ))= S.^FK ^ p) 

where E^ CM is the transverse energy and <fi the azimuthal angle for each EFO, i. The 
value is calculated in bins of rj ncM , semi-inclusively, and for different minimum cuts on 
the transverse energy of the EFO, iffi ™ . Equation @ is also used for the determination 
of the moments in theoretical calculations, in which the sum is either over final-state 
hadrons or partons. The mean values are not expected to be sensitive to uncertainties 
in parton distribution functions, fragmentation functions, and calorimeter energy scale, 
since such effects contribute to both the numerator and the denominator. 



4 Correction procedure 

Monte Carlo (MC) events were used to correct the data for detector acceptances. For all 
generated events, the ZEUS detector response was simulated in detail using a program 
based on Geant 3.13 [20]. 

Neutral current events with electroweak radiative corrections were simulated with the 
Lepto 6.5.1 [21] program interfaced to Heracles 4.6.1 [22,23] via the Djangoh 1.1 
program [24,25]. High-order QCD processes were simulated using the MEPS option of 
Lepto. A second sample of MC events was generated with Ariadne 4.12 [26], where 
the QCD cascade is simulated with the colour-dipole model. In all cases, the events were 
generated using the CTEQ4D parton density parametrisation [27] of the proton. The 
final-state partonic system was hadronised using the Lund string model as implemented 
in Jetset 7.4 [28]. 

For a given bin, j, in r] RCM and F(n0 HCM ), two correction factors for F(n0 HCM ) were 
derived from MC as the ratio of values determined using hadrons and EFOs: 
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[l^i^T,i * {n<Pi jj hadij [l^i^TA J had>j 



5 " [Ei E«™F (^HCM }] ^ ' ' - ^ ^HCM] MC ^ ' 

where Cj corrects the trigonometric function weighted with and Dj corrects the sum 
of Et in the bin, j. The corrected data for a bin in rj M is then given by: 



(F (W> HCM )> 



The overall correction factors are about 10% and arise mainly from undetected hadrons. 

For this approach to be valid, the uncorrected energy flow in the data must be well 
described by the MC simulations at the detector level. This condition was satisfied [29] 
by both the Lepto and Ariadne simulations in the rj HCM and E^ CM regions under 
investigation. The samples of Lepto events were used for the final corrections. 

Systematic uncertainties of the azimuthal asymmetry were determined by varying the 
event selection cuts within their reconstruction resolution and the total systematic uncer- 
tainty was taken as the sum in quadrature of the individual systematics. The dominant 
contributions originated from the following sources (maximal deviations for (cos</>) are 
shown in parentheses): 

• the use of the Ariadne MC model to correct the data (0.017); 

• varying the cut on the pr of the final-state objects from 0.15 GeV to 0.2 GeV to 
estimate the effect due to low-px tracks (0.013); 

• the inclusion of energy deposits consistent with being a photon for the calculation of the 
hadronic angle, used in the double angle method [16], and therefore the transformation 
to the HCM frame (0.009). 

The effect of the variation of other selection cuts was negligible. 



5 QCD calculations 

The data were compared to the MC programs Lepto and Ariadne, described in the 
previous section, and to a next-to-leading-order (NLO) prediction. In these two MC 
programs, the azimuthal-angular distribution is implemented according to the first-order 
QCD matrix elements for QCDC and BGF. The parton-shower and soft-matrix-element 
events are flat in azimuth. The NLO predictions were calculated using the dipole factori- 
sation formulae [30] implemented in the Disent program [30,31]. The calculations used a 
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generalised version of the subtraction method [32] and were performed in the massless MS 
renormalisation and factorisation schemes. The azimuthal-angle distribution is calculated 
to NLO by the Disent program through the use of the photon leptonic current in the 
amplitudes for the QCDC, BGF and 0(a 2 ) processes. This program contains neither Z° 
exchange nor hadronisation effects. The following settings were used as defaults for Dis- 
ent: the number of flavours was five, the factorisation and renormalisation scales were 
Hf,r = Q, and the parton distribution function was CTEQ3M [33]. Samples of events 
from Lepto 6.5.1 were used to correct the NLO QCD calculations for Z°-exchange ef- 
fects, hadronisation and undetected particles due to the requirements on and 9 of the 
produced hadrons. 

The uncertainty in the Disent predictions was estimated by changing the following: 

• the renormalisation and factorisation scales were individually changed to Hf,r — Q/% 
and 2Q; 

• the parton distribution functions CTEQ4M [27], CTEQ5M [34], were used; 

• the correction for Z°-exchange effects, hadronisation and undetected particles was 
repeated with the Ariadne 4.12 MC. 

The above uncertainties were at most 0.004 in both (cos0) and (cos 20). They were 
added in quadrature and are displayed in the figures as shaded bands around the central 
prediction. The sensitivity to different gluon distributions was checked by using the 
MRST99 [35] parton distributions functions with an increased or decreased gluon density; 
the differences were negligible. 



6 Results 

Azimuthal asymmetries have been measured in NC DIS events with the requirements: 
100 < Q 2 < 8000 GeV 2 , 0.2 < y < 0.8 and 0.01 < x < 0. 1 for hadrons with p T > 0. 15 GeV 
and 9 > 8°. The mean values of cos0 HCM , cos20 HCM , sin0 H and sin20 HCM are shown 
in Fig. Eland given in Tables ^ and |2 as a function of r/ HCM . The data are compared with 
predictions from MC models and from NLO QCD as described in the previous section. 

Figure El shows that the value of (cos0 HCM ) is negative for r/ HCM < —2 but becomes 
positive for larger r) nCM . This is in disagreement with both MC predictions, which are less 
negative for r^ HCM < —2 and remain negative for larger ^ HCM . The measured (cos20 HCM ) 
values are consistent with zero for ?^ HCM < —2 but are positive for higher values of ry HCM . 
This is consistent with the expectations from both Lepto and Ariadne. 

The NLO QCD predictions, corrected for hadronisation, agree better with the experi- 
mental values for (cos0 HCM ) than do the MC predictions. The predictions from NLO 
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QCD are more negative for ^ HCM < —2 than those from the MC generators and also have 
positive values for larger ^ HCM . However, the NLO calculation still fails to describe the 
magnitude of the asymmetry in the data. The comparison with NLO QCD was also made 
at higher pr, greater than 1 GeV (not shown). The Monte Carlo was used to correct the 
NLO for this cut; although the final correction for (cos0 HCM ) was small, the correction 
for hadron removal was large. However, the comparison with the data was qualitatively 
the same as when the more inclusive cut was used. The disagreement between data and 
NLO suggests that higher-order calculations may be necessary to describe this distribu- 
tion fully. Inclusion of higher orders through a resummation of large logarithmic terms is 
expected [6] to give an improved description compared to that of LO for —5 < rj < —3. 
However, the description is not significantly better than for the other predictions. For 
(cos20 HCM ), the NLO and MC predictions are similar and describe the data reasonably 
well. 

Figure H3 shows that the values of (sin0 HCM ) and (sin20 HCM ) are small. A deviation of 
(sin HCM ) from zero at the level of three standard deviations is observed. The mean 
values are expected to be at least an order of magnitude smaller than the (cos</> HCM ) 
term [5]. The values of (sin20 HCM ) are consistent with zero. None of the theoretical 
models include predictions for (sin0 HCM ) or (sin20 HCM ). 

To investigate the effect of the minimum transverse energy cut, on the asymme- 

tries, the event sample was subdivided into three regions of r] ncM : —5 < r] ncM < —2.5, 
-2.5 < r] HCM < -1 and -1 < r] HCM < 0. For £™ = 1 GeV, the acceptance is ap- 
proximately 100%. Below this value, some hadrons are removed by the pt > 0.15 GeV 
requirement mainly in the region —2.5 < r] UCM < —1. The data are shown in Fig. 0J and 
given in Tables El and 0] compared to the predictions from Lepto and Ariadne MCs. 
As stated previously, NLO QCD predictions for higher -E^min have large corrections for 
hadron removal. 

The first region —5 < r] ncM < —2.5 is part of the current region in DIS defined in the 
Breit frame as ^ Brcit « ry HCM + 2 < 0; in this region the main contribution to the azimuthal 
asymmetry comes from QCDC and arises from hadrons from quark fragmentation (Fig.|2J). 
This region was investigated in the first ZEUS analysis of azimuthal asymmetries [8] using 
charged hadrons. The data from the current analysis, shown in Fig.HJ confirm, with higher 
experimental precision and to higher E^^^, that the value of (cos0 HCM ) is more negative 
than expected from MC predictions. The (cos20 HCM ) values are small and consistent 
with zero and also in agreement with both Lepto and Ariadne. 

The region —2.5 < r] ncM < — 1 is still dominated by QCDC events but the contribution 
from BGF events increases. The ZEUS analysis of azimuthal asymmetries measured 
using jets [9] was based on hadrons from this region of phase space. In that analysis, 
a large positive value of (cos20 HCM ) was measured, whereas the (cos(/> HCM ) value was 



S 



consistent with zero. The (cos20 HCM ) values agreed with the NLO QCD prediction and 
were inconsistent with the LO prediction. The results presented here in Fig. 0] confirm, 
with higher experimental precision, a small value of (cos0 HCM ) and positive values for 
(cos20 HCM ) for all E^™- The MC predictions of Lepto and Ariadne are in good 
agreement with the data. 

The third region, — 1 < ^ HCM < 0, is populated roughly equally by hadrons from QCDC 
and from BGF processes. This measurement extends the kinematic region over those 
presented previously. The results are presented in Fig. 0J The (cos(/> HCM ) values are 
positive, contrary to MC predictions, whereas the (cos20 HCM ) values are positive and 
in agreement with MC predictions. These trends persist up to the highest E^ ( ^ l values 
measured. 



7 Summary and conclusions 

The azimuthal asymmetries in deep inelastic scattering have been measured in the hadronic 
centre-of-mass frame at HERA for a selected sample of neutral current events with 
100 < Q 2 < 8000 GeV 2 , 0.2 < y < 0.8 and 0.01 < x < 0.1. An energy-flow analysis 
method was used, which permitted the use of both neutral and charged hadrons and 
extends the phase space over previous measurements. 

Azimuthal asymmetries have been investigated as a function of hadron pseudorapidity, 
7/ HCM . The value of (cos0 HCM ) is negative for rj RCM < —2 but becomes positive for 
larger ry HCM . The distribution is not well described by the MC predictions of Ariadne 
and Lepto. Although the predictions from NLO QCD describe the data better than 
do the MCs, they still fail to describe the magnitude of the asymmetries. This suggests 
that higher-order calculations may be necessary to describe these data. However, the 
predicted values of (cos20 HCM ) in both MC models and in NLO QCD agree with the 
data. A deviation of (sin0 HCM ) from zero at the level of three standard deviations is 
observed. The values of (sin 2<j) HCM ) are consistent with zero. 
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Table 1: The values o/(cos0 HCM ) and (cos 2</> HCM ), calculated using the energy- 
flow method as in Eq. as a function of hadron pseudorapidity, 7y HCM . They 
were obtained in the HCM frame for the kinematic region 100 < Q 2 < 8000 GeV 2 , 
0.01 < x < 0.1 and 0.2 < y < 0.8 for hadrons with p T > 0.15 GeV and 9 > 8°. The 
quantities 5 gtat and d~ syst are respectively the statistical and systematic uncertainties. 
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Table 2: The values o/(sin0 HCM ) and (sin20 HCM ), calculated using the energy- 
flow method as in Eq. as a function of hadron pseudorapidity, r] KCM . They 
were obtained in the HCM frame for the kinematic region 100 < Q 2 < 8000 GeV 2 , 
0.01 < x < 0.1 and 0.2 < y < 0.8 for hadrons withp T > 0.15 GeV and 9 > 8°. The 
quantities 8 stat and d~ syst are respectively the statistical and systematic uncertainties. 
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Table 3: The values of (cos0 HCM ) calculated using the energy-flow method 
as in Eq. as a function of hadron minimum transverse energy, E^^f^. 

They were obtained in the HCM frame for the pseudorapidity intervals —5 < 
77 M < — 2.5, —2.5 < jj HCM < — 1 and —1 < r] HCAI <0 in the kinematic region 
100 < Q 2 < 8000 GeV 2 , 0.01 < x < 0.1 and 0.2 < y < 0.8 for hadrons with 
Pt > 0.15 GeV and > 8°. The quantities 5 sta t and S syst are respectively the 
statistical and systematic uncertainties. 
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Table 4: The values of (cos0 HCM ) calculated using the energy-flow method 
as in Eq. as a function of hadron minimum transverse energy, E^^f^. 

They were obtained in the HCM frame for the pseudorapidity intervals —5 < 
rj HCM < - 2.5, -2.5 < rj HCM < - 1 and -1 < r] HCM <0 in the kinematic region 
100 < Q 2 < 8000 GeV 2 , 0.01 < x < 0.1 and 0.2 < y < 0.8 for hadrons with 
Pt > 0.15 GeV and > 8°. The quantities 5 st at one? G- syst are respectively the 
statistical and systematic uncertainties. 
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Figure 1: The definition of the azimuthal angle <p either in the HCM or the 
Breit frame. The incoming electron is denoted by e, the scattered electron by e' , the 
exchanged virtual photon by 7* and the outgoing hadron or parton by h. 
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Figure 2: (a) The fraction of BGF (dashed line), QCDC (full line) and QPM (dot- 
ted line) processes as a function of pseudorapidity, rj HCM ; in the H CM frame for the 
energy-flow method, (b) For the QCD Compton process, the quark and gluon con- 
tributions as a function ofrj . These predictions were taken from Lepto 6.5.1 
and are shown for the kinematic region 100 < Q 2 < 8000 GeV 2 , 0.2 < y < 0.8 and 
0.01 < x < 0.1 for hadrons with p T > 0.15 GeV and 9 > 8°. 
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Figure 3: The values of (cos (j> HCM ) , (cos20 HCM ), (sin^ HCM ) and (sin20 HCM ) ; 
calculated using the energy-flow method as in Eq. as a function of hadron 
pseudorapidity, r] HCM . They were obtained in the HCM frame for the kinematic 
region 100 < Q 2 < 8000 GeV 2 , 0.01 < x < 0.1 and 0.2 < y < 0.8 for hadrons with 
Pt > 0.15 GeV and > 8°. The inner error bars are statistical uncertainties, the 
outer are statistical and systematic uncertainties added in quadrature. The NLO 
QCD predictions of Disent (solid line), with its associated uncertainty (shaded 
band), corrected for hadronisation and hadron losses (see text), the predictions of 
Lepto 6.5.1 (dotted line), and the predictions of Ariadne J^.12 (dashed line) are 
shown. 
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Figure 4: TTie values of (cos(j) HCM ) and (cos2cj) HCM ) , calculated using the energy- 
flow method as in Eq. ffi), as a function of hadron minimum transverse energy, 
^Tmin- They were obtained in the HCM frame for the pseudorapidity intervals 
-5 < n HCM < - 2.5, -2.5 < n HCM < - 1 and -1 < <q HCM <0 in the kinematic 
region 100 < Q 2 < 8000 GeV 2 , 0.01 < x < 0.1 and 0.2 < y < 0.8 for hadrons 
with pt > 0.15 GeV and 9 > 8°. The inner error bars are statistical uncertainties, 
the outer are statistical and systematic uncertainties added in quadrature. The 
predictions of Lepto 6.5.1 (solid line) and of Ariadne (dashed line) are 
shown. 
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